Creating vortex rings and three-dimensional skyrmions in Bose-Einstein condensates 
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We propose a method of generating a vortex ring in a Bose- 
Einstein condensate by means of electromagnetically-induced 
atomic transitions. Some remnant population of atoms in a 
second internal state remains within the toroidal trap formed 
by the mean field repulsion of the vortex ring. This population 
can be removed, or it can be made to flow around the torus 
(i.e. within the vortex ring). If this flow has unit topologi- 
cal winding number, the entire structure formed by the two 
condensates is an example of a three-dimensional skyrmion 
texture. 

03.75.Fi,05.30.Jp 

With the recent controlled creation of vortices in 
trapped atomic Bose-Einstein condensates (BECs), and 
evidence for vortex nucleation above a critical flow veloc- 
ity 1^, the dynamics of topological structures in weakly 
interacting superfluids has become an active subject of 
experimental research. In this paper we propose meth- 
ods of preparing a vortex ring within a BEC, allowing 
the controlled study of a vortex structure away from sur- 
face effects Q. The vortex ring is created by means of 
driving electromagnetic (em) fields inducing transitions 
between internal atomic levels. The em field couples to 
the relative phases between the levels, in such a way as 
to transfer angular momentum to the atoms. With an 
appropriately chosen superposition of em fields the lo- 
cal axes of this effective rotation form a closed circular 
loop, producing a vortex ring in the BEC. We also con- 
sider the stability of the vortex ring, and identify a way 
to use a second atomic component of the BEC to in- 
crease the stability of the ring against collapsing under 
its own string tension. Conversely, the vortex ring acts as 
a trap for the second component, providing purely atom- 
optical confinement of cold atoms. Finally, the combined 
two-component structure may be identified as a three- 
dimensional (3D) skyrmion Q, shown in Fig. ^ This 
nonsingular texture is of interest in its own right, as a 
topological structure beyond the simple vortex. 

Topological defects such as vortices and monopoles 
possess cores at which the order parameter is singular, 
and indeed such defects are characterized by winding 
numbers defined on paths or surfaces that enclose the 
singular cores. The simplest topological defect relevant 
to superfluid physics is that of a current flowing around 
a closed ID path parametrized by the angular coordi- 
nate (j). The simple mapping a{(j)) — cj), from physical 
space into order parameter space, where <j) is defined on 



a ID circle, has 'winding number' one, which will be un- 
changed if a((/)) is continuously deformed. Topologically 
nontrivial configurations are also possible without singu- 
lar cores: Textures are defined by the way in which the 
compact order parameter space is 'stretched over' physi- 
cal space. An SU(2) order parameter, for instance, takes 
values on a sphere. One can 'puncture' this sphere, 
stretch the resulting hole to infinity, and so 'spread' the 
order parameter space over a 2D plane. Such a mapping, 
or any continuous deformation of it, is known as a 2D 
(or 'baby') skyrmion. Interpreting the order parameter 
sphere as the spin orientation of a spinor condensate, one 
can realize a 2D skyrmion as a coreless point vortex (or in 
3D, a coreless vortex line) [^[|^, such as has been created 
experimentally at JILA m. 




FIG. 1. Schematic illustration of a 3D skyrmion (left), as 
a vortex ring containing a superflow, and a constant surface 
plotting of the 3D skyrmion (right) obtained from the nu- 
merical simulations. We show the core region of the vor- 
tex line and the ring by displaying the atom densities for 

IV'i|' = IV'2|' = i.2xio-V;3. 

The full, 3D skyrmion results from taking an order 
parameter space and stretching it over in a similar 
way [||. An order parameter space is afforded most 
simply by a two-component BEC, whose interactions ef- 
fectively fix the total density IV"! P + l''/'2p = P of the two 
complex macroscopic wave functions iI)It^2- The low en- 
ergy degrees of freedom therefore reside on a 3-sphere 
of angles < a, /3 < tt, < 7 < 27r (a 3-sphere is a 
3D subset in 4D, such that Ylti=i •'^n — ^^)- For spher- 
ical polar coordinates in physical space, a 3D skyrmion 
is given by any continuous deformation of the mapping 
■y{r,9,(f>) — <j), (3{r,9,(l)) — d,a{r,9,(p) = \{r) for mono- 
tonic A with A(0) — 0; A(oo) — it. The topologically 
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invariant winding number 

W = J (fx sin^ a sin (3 det (^^^ , (1) 

which obviously vanishes for the ground state of the two- 
species BEC, equals one for the skyrmion. (The topolog- 
ical invariance of W is not hard to understand: because 
the determinant in the integrand is precisely the Jaco- 
bian for the change of variables from = {x, y, z) to 
a' — (a,/3,7), the integral is simply the volume of the 
unit 3-sphere, regardless of the precise form of a*(a;^).) 

We can obtain a physical picture of the 3D skyrmion 
by writing out the macroscopic wave functions: 

/■0i(r)\_ rrpr f -isinasmf3exp{ij)\ 
I -02(1") y V ) \ cos a — i sin a cos /3 J 

-.r^( -isin[A(r)]sin6'exp(i(/)) \ , . 

-ypy^iycos[\{r)]-isin[x{r)]cose) ■ 

With this realization of the skyrmion, ji/'iP vanishes on 
the z axis and at infinity, and is concentrated in a toroidal 
region. On the other hand, |'02p vanishes on the circle 
9 = ^,1" = A~^(^). Expanding 'ip2 about any point on 
this circle, we find V'2 oc Sr+i69, indicating that the nodal 
circle of ■02 is a vortex line. Hence component -02 forms 
a vortex ring, within whose core ipi resides and flows az- 
imuthally with winding number one. Of course, any con- 
tinuous deformations of Eq. (^) will still be a skyrmion, 
and the two-fluid interpretation of the skyrmion will al- 
ways provide a useful intuitive picture. Moreover, unlike 
in the case of wave function monopoles |^ , changing the 
atomic basis that defines the two fluids does not lead to 
a different picture of the skyrmion, but simply produces 
the same combination of vortex ring and current, with a 
different spatial orientation. 

We propose a method of creating a 3D skyrmion in a 
trapped BEC using em fields, in a scheme motivated by 
viewing the skyrmion as this combination of vortex struc- 
tures. Since several good methods have been proposed 
for generating simple vortices, and some of these have 
already been successfully implemented, engineering the 
state of 01 presents no new difficulties of principle. Cre- 
ating the vortex ring state in 02 is more of a challenge. 
Although vortex rings may be formed through instabili- 
ties, as of a dark soliton in an insufficiently narrow trap 
Q , a more controlled method of preparing them may also 
be desirable, in order to study their rich range of motion 
and excitation. 

We consider em-induced coherent scattering of atoms 
between different internal levels. By superposing differ- 
ent driving em fields we can construct a Rabi amplitude 
with a spatially-dependent phase profile. In particular, 
nodes of the em field amplitude may be topological sin- 
gularities of the em phase. The coupling of the em field 
to the phase of the coherent matter wave allows the topo- 
logical singularities of the em fields to be imprinted on 



the matter field. We demonstrate the preparation pro- 
cess of topological objects by numerically studying this 
coherent matter- wave dynamics. 

The dynamics of the BECs with the Rabi coupling 
between the levels \i) and \j) follows from the coupled 
Gross-Pitaevskii equation (GPE) 

ihljj, = iH° +6,+J2 ^^k\^k\^)^^ + ^^^0 ■ (3) 
k 

Here the kinetic energy and the trapping potential are 
introduced in H^: 

^ -|1 + Imu^Ux' + ajy' + (3fz') . (4) 

We have also defined the interaction coefficients Kij = 
iTTh^aijN/m. Here an denotes the intraspecies scattering 
length in internal level \i) and Uij (i ^ j) stands for the 
interspecies scattering length. The Rabi frequency, 0(r), 
describes the strength of the coupling between the two 
internal levels. The total number of BEC atoms, the 
atomic mass, and the detuning of the cm fields from the 
resonance are denoted by N, m, and 5j, respectively. 

Several papers [p^j-pslj?! have proposed methods of in- 
jecting vorticity into BECs by means of the interaction 
between em and matter fields. In our adaptation of these 
methods to create a 3D skyrmion, the simplest part is the 
step of producing the 'inside' component -01 with a quan- 
tized circulation about the z axis. This can be achieved 
by transferring population from an initial internal state 
0) using the Rabi frequency noi(r), 

f2oi(r) ~ f7o[sin(fcix) — isin(fci?/)] cos(fc2z) , (5) 

with a similar set of parameters as in Ref. p5| ]. In par- 
ticular, fcf ^ and fc^^ are both larger than the size of the 
sample. The purpose of the cos(fc2z) is to confine 0i 
along the z direction. Note that in this scheme the Rabi 
frequency could represent a two-photon transition, with 
one photon from two orthogonal standing waves in the 
xy plane and the other from a standing wave along the z 
direction. Alternatively, the cos(fc2z) could simply rep- 
resent collimation of the xy beams, which would drive a 
single photon transition by themselves. The topological 
singularity in the em field at x = y = results in a single 
quantum of circulation about this axis, in the component 
ipi, as displayed in Fig. ||. One may also generate addi- 
tional terms in Eq. (|5|). For instance, the standing wave 
along the z axis could be further coUimated in the xy 
plane, limiting the spatial extent of "01- 

We propose two different methods of creating the more 
complicated part of the skyrmion, the closed vortex ring 
in 02. The first is an adaptation of the method of Ref. [Q, 
and makes use of two internal levels and the circulating 
state ipi once it has been generated. To add the vortex 
ring requires a Rabi field coupling from state |1) to state 
2), which is a phase-coherent superposition of a standing 
em field and a Gaussian beam, both along the z axis: 



2 



Quir) = r!o [1 - - exp [-pV^'l + lAsinihz)] , (6) 

with A = 2r]/{k3() and p = {x'^ + y^Y^^. We are inter- 
ested in kslzl ^ 1, p <C ^, and r/ 1 in such a way that 
fc3(p — 77^)1 <C 1. Then we obtain 

ni2ir)/no ^ -V^ +^Ak3Z 

^Aks[ip-vO+^^]- (7) 

The phase of rii2(r) has a form of the quantized circula- 
tion with unit winding number around the closed ring p — 
rjS^ and z = 0. Alternatively, the Gaussian beam resulting 
in the exponential term in Eq. (^) may also be replaced 
by two standing em fields: [cos(2x/^) + cos{2y /£_)]/ 2 . 

For this first method, the field of Eq. is tuned away 
from resonance. A localized shifting field is then applied 
as m Rcf. (1), to 'turn on' the transition in a small region, 
which is then swept around the trap to create a full vor- 
tex ring, and also undo the azimuthal phase winding of 
the |1) atoms, so that the population shifted into |2} will 
have no angular momentum about the z axis. While this 
scheme takes advantage of an already demonstrated tech- 
nology, it has the disadvantage that it tends to produce 
a highly excited skyrmion state, which may (as discussed 
below) then evolve self-destructively. 




FIG. 2. A 3D skyrmion in the xy plane. We show the 
density \tlj{x,y)\^ and the phase — aig[tlj{x, y)] profiles of the 
wave function in the xy plane. Here ip{x,y) = J dzil){x,y, z). 
The density profile (upper left) and phase profile (lower left) 
of atoms in level |1) display a single vortex line with 2-k phase 
winding around the core, and vanishing density at the center 
of the core. The vortex line occupies a toroidal region where 
the density of the vortex ring of atoms in level |2) (upper 
right) vanishes. The phase distribution of the vortex ring 
(lower right) in the xy plane is fiat. 

We therefore propose, in our second method, to couple 
atoms with em fields from level |0) to level |2), at the 
same time as we generate the simple vortex state in level 
1). We need a Rabi field coupling from |0) to |2), of 
the same spatial form as of Eq. (H), as well as the 
fioi field of Eq. (||). This scheme requires control of an 
additional internal state |0), but with a weak coupling it 



allows adiabatic transfer of the atomic populations [ [l2[ , 
producing a ground state skyrmion. 

For computational simplicity we demonstrate the 
three-level preparation process by considering two short 
nonadiabatic Rabi pulses: The first Rabi pulse of the 
form (^) generates the vortex ring in level |2) from the 
initial ground state population in level JO), and then the 
second pulse with the configuration (|g) generates the 
trapped current in level |1) from the remaining atoms in 
level |0). Our numerical approach uses three states and 
can still generate an excited skyrmion, but it is amenable 
to fully 3D simulation, and our simulation indicates that 
the method can work. A careful study of an adiabatic 
population transfer will require a realistic model for the 
dissipation of the BECs, as well significant computational 
resources, and must be left for future investigation. 

In our numerical simulations we assume that the traps 
are isotropic, = /3i = 1, and that the trapping fre- 
quences are equal for all components. We choose the 
wave number of the em fields ki — k2 — k^ = 2tt/{81), 
with I = [?i/(mw)]i/2, S, = 0, r] ^ 0.04, ^ = 30/, and 
Kij — lOOhujl^. The number of atoms is determined by 
the previous relations as = 251 /{tto). 




FIG. 3. A 3D skyrmion in the yz plane. We display the 
density \tpix = 0,y,z)f and the phase = arg[?/)(a; — 0,y,z)] 
profiles of the wave function. The density of the single vor- 
tex line in level |1) (upper left) occupies the toroidal region 
surrounded by the vortex ring in level |2) (upper right). The 
phase profile of the vortex ring (lower right) displays the quan- 
tized rotation around the ring. 

We integrate the GPE on a spatial grid of 128"^ points, 
using a split-operator method. The stability of the non- 
linear evolution is improved by simultaneously propagat- 
ing multiple copies of the wave function. Relatively small 
number of time steps (less than one thousand) provide 
a sufficient accuracy making the integration executable. 
Even with up-to-date workstations the optimized simula- 
tions can take several hours. Figures represent the 3D 
skyrmion in the case where we have first created the vor- 
tex ring in level |2) by applying the Rabi frequency rii2(r) 
[Eq. with rio = 5000w for < tcu < 0.04. This is 
followed by creating a vortex line in level |1) by the driv- 
ing f}oi(i-) [Eq. (|)], with Qo = 80uj for 0.04 < tuj < 0.05. 
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With the chosen set of parameters a large population re- 
mains in [0). This remnant population may be coupled 
out of the trap immediately after the em pulses as in the 
recent vortex experiments Due to the shortness of 
the pulses the interactions of the atoms in level |0) do 
not have time to affect the structure of the skyrmion. 
On the other hand, for an efficient adiabatic population 
transfer the nonlinearity can be a problem. However, the 
techniques to overcome this were developed in Ref. . 

The effort of creating such an intricate structure as a 
vortex ring, let alone a 3D skyrmion, is obviously only 
justified if its structure can then be observed. As with the 
filled vortex lines of JILA, the interior BEG in state |1} 
will enable detection of the density and phase profile of 
the skyrmion. An additional advantage of the skyrmion, 
as opposed to an empty vortex ring or to a vortex ring 
filled with a nonrotating interior BEG, is that it is less 
unstable against shrinking to zero radius. The kinetic 
and interaction energy of a vortex line gives it an ef- 
fective string tension tending to shorten the vortex line, 
and in the case of a vortex ring, the ring's radius R (not 
its core thickness!) is energetically unstable to shrink- 
ing. In an infinite, homogeneous background, and in 
the absence of dissipation, the first-order dynamics of a 
vortex line means that a vortex ring will not actually 
shrink, but will instead move along its axis, at a velocity 
determined by its radius. But if its radius approaches 
its core thickness, a vortex ring can be annihilated (it is 
not topologically stable) , and an excited vortex ring may 
also possibly twist or pinch, and break up into smaller 
rings. Moreover, motion at a constant speed cannot con- 
tinue in a finite, inhomogeneous BEG. Dissipation is also 
present in experimental conditions, and although straight 
vortices have been found to have a long lifetime before 
floating out of the BEG cloud, the energy slope due to 
string tension may be much steeper than that due to 
vortex buoyancy (on the healing length scale rather than 
trap scale). Filling the vortex core lowers the string ten- 
sion, and hence extends the lifetime of the filled vortex 
structure in comparison with that of an empty ring. 

Although in other field theories skyrmions are often 
unstable against shrinking to zero size without the topo- 
logical change of defect nucleation, the separate conser- 
vation of both atomic species in our case forbids this de- 
cay channel. This does not prevent the vortex ring from 
collapsing, but it does mean that the component |1) fill- 
ing the ring cannot be eliminated. Furthermore, since 
the filling has one unit of angular momentum, there is a 
centrifugal barrier inhibiting the shrinking. As it tends 
to pin the vortex ring down, the ring may only collapse 
via the tunneling of atoms through ipi. Although the 
vortex rings are dynamically stable [|l6| , numerical evolu- 
tion in imaginary time shows that filled vortex rings, and 
even skyrmions, are still energetically unstable, on long 
enough (imaginary) time scales. Skyrmions are topologi- 
cally stable as long as the order parameter is well defined 
everywhere; but they can be destroyed by the nucleation 
of a defect. Our simulations indicate that, while the en- 



ergy decreases slowly but monotonically, the vortex ring 
in 1^2 can pull through the state |1) if its inner radius 
approaches the healing length, so that a healing-length 
radius circle of zero total density nucleates. However, the 
decay via the tunneling of atoms through -01 becomes in- 
creasingly slow with the larger occupation in |1) jl^. A 
vortex ring with a nonrotating filling can, if its radius 
becomes too small, annihilate itself without significantly 
lowering the total density at any point. 

Without dissipation, numerical evolution in real time 
shows that skyrmions as well as empty vortex rings do 
move through the background cloud, but instead of sim- 
ply exiting the cloud as they approach its surface, they 
generate a background flow in the rest of the cloud which 
is able to draw them back towards the center. Much more 
extensive numerical studies will be needed to determine 
the stability of the resulting oscillation. Further analyt- 
ical and numerical studies of skyrmion and vortex ring 
dynamics are clearly warranted, but in this paper we have 
shown that engineering and study of these advanced BEG 
structures should be feasible with current techniques. 
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